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Studies  of  the  dynamics  of  several  brain  enzymes  and  membrane  receptors  indicate  that  the  pres¬ 
ence  of  hierarchical  modes  in  power  law  spectra  is  associated  with  the  capacity  for  signal-sensitive  con¬ 
formational  transitions.  When  transient,  this  transition  to  differentiability  is  characteristic  of  the  system 
in  action;  when  long-lasting,  the  mode  locked  periodicity  leads  to  receptor  protein  denaturation  and 
desensitized  preparations.1-3  Multiplicity  in  substrate-velocity  curves  (iterative  saturation  plateaus)  of 
allosteric  proteins  demonstrated  comparable  multimodality  in  velocity,  V,  across  increasing  substrate 
concentrations,  S,  in  a  staircase  of  half-maximal  rate  constants  in  addition  to  adaptive  renormalization 
as  many  smaller  steps  became  fewer  larger  ones  when  the  concentrations  of  regulatory  ligands,  r,  were 
varied  (Fig.  I).3-5 
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Fig.  1  Ligand-induced  renormalization  of  iterative  substrate  saturation  plateaus  in  catalytic  velocity  of  allosteric 
proteins. 

We  note  the  loss  of  multimodal  complexity  in  enzymatic  and  receptor  function  away  from  a  qua- 
dratically  maximal  value  of  r.  The  problem  involves  understanding  both  structural  stability  with  respect 
to  the  underlying  equations  (not  initial  values)  and  phase  transitions  in  these  dynamical  hierarchies. 
Evolutionary  selection  of  restricted  ranges  of  biological  parameters  removes  genericity,  density,  or  even 
positive  measure  as  necessary  to  such  explorations.  For  example,  our  theory  says  that  the  signal- 
sensitive  allc'teric  protein  lives  in  the  e -neighborhood  of  some  complex  boundary  on  the  hyperbolic 
plane.  In  real  parameter  space,  these  are  the  irregular  border  regions  of  competing  attractors  demon¬ 
strating  homoclinic  tangencies  and  self-similar  scaling. 

A  Thermodynamic  Law  for  Signal-Sensitive  Proteins 

The  universality  of  linear  free-energy  relationships  between  250  and  315°K  of  both  rate  and 
equilibrium  constants  for  a  wide  range  of  water-solvated  organic  chemical  reactions,6  including  enzy¬ 
matic  catalysis  and  protein  ligand  binding7  by  biological  macromolecules  is  seen  in  the  form  of  a  nor¬ 
malizing  factor  called  compensation  temperature,  a,  which  rescales  the  temperature  at  which  a  class  of 
reactions  are  run,  leaving  invariant  the  slopes  of  the  linear  Gibbs  free  energies  of  interaction,  a  has 
also  been  called  the  isokinetic  or  iso-equilibrium  temperature  because  Arrhenius  or  Van't  Hoff  plots 
intersect  at  1/a  when  temperature  compensation  is  present.  The  molecular  mechanism  underlying 


those  findings,  observed  most  systematically  in  studies  of  the  effect  of  aliphatic  substitution  on  ioniza¬ 
tion  rates  of  organic  acids  using  the  Hammett  equation,  has  been  called  the  hydrophobic  effect.8  Water 
adjusts  its  geometry  to  maximize  the  number  of  intact  hydrogen  bonds  with  itself  to  neutralize  the 
unbonded  electron  pair  of  oxygen  by  forming  clathrate-like  convex  cages  around  flickering  aggregates  of 
self-seeking,  non-soluble  organic  moieties  while  in  a  dynamical  struggle  to  minimize  lost  solvent 
entropy.  The  hydrophobic  effect,  a  function  of  the  aqueous  surface  (•—  l2),  scales  the  mass  (—  I3)  of 
coherent  aggregates  of  polypeptide  polymers  as  mV3.  This  estimate  has  been  confirmed  in  more  than 
30  proteins  up  to  50,000  Daltons.9  The  hydrophobic  effect  as  a  hyperbolic  interfacial  free  energy  scales 
self-similarly  across  the  hierarchies  of  peptide  clusters  in  water  reflected  in  the  fractional  exponential 
renormalization  of  their  temperature-dependent  kinetics.  In  addition,  the  isomorphism  of  their  scaling 
exponents  in  Euclidean  space  R2  and  complex  space  C,  x-ray  crystallographic  and  density-of-states  spec¬ 
tra  respectively  (see  below),  is  consistent  with  a  thermodynamic  rather  than  a  physical  anatomy  of  these 
macromolecules. 

Water’s  unusual  propensity  for  dynamical  self-rearrangement  into  states  that  are  both  attractive 
and  repelling  eliminates  negative  values  for  the  Gibbs  free  energy  of  conventional  equilibria.  Decom¬ 
position  of  the  flow  into  stable  and  unstable  manifolds  makes  rigorous  thermodynamic  treatment  inap¬ 
plicable.  That  AG  >  0  in  aqueous  solutions  of  small  as  well  as  polymeric  organic  molecules  can  be 
called  enthalpy-entropy  compensation;  AH  =■  a  AS  -f  A  Ga  implies  that  in  the  neighborhood  of  tem¬ 
perature  T°a,  AH  is  paralleled  by  a -renormalized  AS,  leaving  the  free  energy  change  AG  negligible  and 
almost  independent  of  the  temperature.  Using  the  second-law  relationship,  AG  -  AH  -  T°AS,  we  can 
rewrite  the  expression  of  enthalpy-entropy  compensation  to  demonstrate  the  role  of  factor  a  with 
respect  to  the  linear  free  energy  as  A  Gy*  *  AH(1  —  T°/a)  +  T°AGa/a  — -  0.  Using  estimates  of 
interfacial  Free  energies  from  surface  area  and  microcalorimetry,  calculations  of  the  countervailing 
forces  of  water  hydrogen  bond  strain  minimization  by  protein  hydrophobic  packing  and  equivalent  loss 
of  polypeptide  chain  entropy  in,  for  example,  hen  egg  white  lysozyme  (AG  —  350  kcal  (AH)  -  350 
kcal  CPAS)  *  0)  exemplify  the  way  protein  dynamics  mirror  the  hyperbolic  stability  of  aqueous  sol¬ 
vent.9  Thermodynamic  intuition  dictating  the  axiom  of  choice  based  on  the  variational  principle  of 
energy  minimization  is  replaced  by  one  of  geometric  limits  on  the  maximization  of  hyperbolic  pressure 
on  a  quadratic  real  or  complex  manifold.10  Allosteric  proteins  rest  on  separatrices  between  basins  of 
attraction,  hierarchically  scaled,  competing  states  of  hydrophobic  aggregation  and  maximum  possible 
chain  entropy  with  well  documented  correlation  times  stretching  from  10” 12  to  103*  seconds.3  One  can, 
for  example,  visualize  the  hierarchically  self-similar,  homoclinic  orbits  and/or  horseshoes  that  consti¬ 
tute  the  strange  point  sets  of  Poincare,  Birkhoff,  and  Smale.  The  transient  global  C°—  C>1  transition 
to  smoothness  represents  the  making  of  functional  connections  between  elements  of  the  hierarchical 
machinery.  This  C°  resting,  sensitive  condition  stabilizes  against  the  spontaneous  emergence  of  C*1, 
coherent  biological  actions  without  the  appropriate  signal.  OnC  this  boundary  involves  the  properties 
of  univalence  and  discreteness.  Signal  transduction  involves  phase  transitions  to  differentiable  dynam¬ 
ics  followed  by  relaxation  to  a  previous  state  of  discretely  mappable  scaling  hierarchies.  Quantitative 
universality  in  macromolecular  dynamics  arises  at  the  boundary  of  differentiability  on  R  and  discrete¬ 
ness  on  C  as  a  "fixed  point."  The  region  of  transition  can  be  defined  by  a  unique  condition,  the  pres¬ 
ence  of  iterative  roots  of  a  Fibonacci  recursion  relation,  on  the  embedding  of  discrete  semigroups  of 
nonlinear  transformations  into  continuous  ones,  P  (■)  —  F(-,  t)  :R  —  R  and  at  a  comparable  boundary 
defined  by  the  transformational  isometries  of  Lobachevskii  space,  T(az  +  b/cz  +  d): L—  L. 

A  continuity  in  the  thermodynamic  state  space  between  the  two  phenomenologically  different 
kinds  of  water,  hydrophobically  structured  and  free,  creates  the  characteristic  family  of  coexistence 
curves  found  in  the  neighborhood  of  second-order  phase  transitions.  In  real  space,  solvent  molecules 
are  interspersed  intimately  with  the  amino  acids  inside  globular  macromolecules  at  a  ratio  of  about  3 
H20:1).  The  lazy  cubic,  "sigmoicT  kinetic  functions  of  brain  enzymes  and  membrane  receptor  prepara¬ 
tions  across  changing  concentrations  of  regulatory  ligands11  exemplify  the  continuous  phase  transitions 
of  allosteric  proteins. 

J.  Willard  Gibbs12  geometrized  the  thermodynamic  state  space  of  a  simple  fluid  in  equilibrium  as  a 
once  differentiable  manifold  with  the  shape  of  a  convex  body  and  a  boundary  upon  which  two  pure 
phases,  such  as  hydrophobically  structured  and  free  water,  meet  at  a  critical  point.  His  phase  rule 
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goveming  components  /,  coexisting  phases  p,  and  dimensionality  d  of  the  manifold  M  of  their  coex¬ 
istence,  i  +  2  —  p  ”  'Ad  >  0,  predicts  that  the  pure  phases  of  the  protein-water  system  will  meet  on  a 
line.  A  cross-section  tnrough  this  entropy  surface  creates  a  two-phase  coexistence  manifold.  At  critical 
values  of  the  parameters,  a  hierarchy  of  metastable  states  emerges  in  the  neighborhood  of  the  critical 
point,  here  represented  by  two  descriptions  for  the  same  set,  superstructured  free  and  superfree  struc¬ 
tured  water  as  extensions  of  convex  coexistence  surfaces  (Fig.  2).  A  Gibbs-like  plot  of  a  cross-section 
through  the  concave  entropy  surface  diagrams  how  correlation  length  diverges  with  proximity  to  the 
critical  point  in  continuous  phase  transitions,13  creating  a  geometrically  proportional,  dilation  symmetric 
hierarchy  with  universal  scaling  exponents  of  the  sort  observed  in  the  relaxation  spectra  of  proteins.3 

Critical  point 
mixture  of  structured 
and  unstructured  water 

Metastable 
superstructured 
free  water 


Fig.  2  Correlation  lengths  along  the  hierarchies  of  convex  co-existence  surfaces  diverge  near  the  critical  point. 

The  intrinsic  geometric  structure  of  R 2  in  interaction  with  a  system  having  cycles  of  many  orders 
creates  such  orbital  complexity  that  its  description  often  requires  Markoff  partition,  symbolic  dynamics, 
and  invariant  measures  such  as  topological  or  measure  theoretic  entropies.14  Physical  (geometric)  intui¬ 
tion  is  lost  as  the  price  for  quantitative  characterization.  Although  the  complicated  behavior  is  gen¬ 
erated  by  deterministic  equations,  its  description  in  entropies  leads  to  even  less  insight  about  the 
mechanism  than  the  assumptions  of  a  describably  constrained  but  purely  random  process  such  as  a 
self-avoiding  walk.  The  focus  on  "sensitivity  to  initial  conditions"  and  "rate  of  information  los£'  has 
generated  nihilism  about  a  physical  understanding  of  the  kind  of  hierarchical  dynamics  that  characterize 
the  behavior  of  signal-sensitive  proteins.1  Perhaps  consideration  of  manifolds  with  other  intrinsic 
geometries  will  be  more  rewarding.  For  example,  a  global  property  of  aqueous  solvent  such  as  its 
hydrophobic,  hyperbolic  pressure  can  be  represented  as  a  general  characteristic  of  the  thermodynamic 
manifold,  its  curvature. 


A  Natural  Metric  for  the  Critical  Dynamics  of  Allosteric  Proteins 

The  intrinsic  geometry  of  the  manifolds  of  thermodynamic  minima  manifests  positive  constant 
Gaussian  curvature,  K  >  0,  a  recurrent  orbit  has  smaller  circumference,  c,  than  it  would  have  on  the 
Euclidean  plane  R2  (Fig.  3a).  In  the  neighborhood  of  a  stable  manifold  of  an  orbit  with  expanding  as 
well  as  contracting  eigenmotions,  the  hyperbolic  Lobachevskii  space  L  with  constant  negative  curvature, 
K  <  0,  stretches  c  beyond  its  length  on  R 2  so  that  for  small  r  its  dependence  on  radius  r  deviates  from 
27 rr  like  c(r)  *■  2-rrr  -  n/3  •  K r3  +  «r3  in  which  «  — *  0  as  r  — *  0  (Fig.  3b). 


a)  b) 


Fig.  3  Manifolds  of  positive  (a)  and  negative  (b)  Gaussian  curvature. 
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This  phase  space  (Fig.  3b)  can  be  represented  by  Lobachevskii’s  geometry  defined  by  replacement 
of  Euclid’s  fifth  postulate  about  parallel  lines  as  seen  in  a  construction  (Fig.  4)  where  two  lines  (here 
symmetric)  through  points  P,F  lie  parallel  to  the  given  line  A,A',  separated  by  their  angle  of  parallel¬ 
ism  <f>  —  f (r)  lying  between  the  limit  lines  of  those  contained  in  <f>  that  converge  to  A, A'  at  infinity.  In 
L  space  the  "equaT  distances  E,E'  and  E',E"  between  parallel  lines  A,A':B,B\  C ,C  appear  to  shrink  at  a 
geometrically  proportional  rate  dependent  on|r  —  r£,  r0  being  in  an  R 2  neighborhood.  A  renormaliza¬ 
tion  of  this  length  that  is  appearing  to  shrink  in  geometric  proportion  along  parameter  r  in  R2  generates 
metric  universality  in  a  geometric  series  that  is  uniquely  both  multiplicative  and  additive,  a  Fibonacci 
recursion  relation,  <Pn  — *  4>',+  l...  andd>"-1  +  4>"  * 


Fig.  4  Lobachevskii’s  geometry  defined  by  the  replacement  of  Euclid’s  fifth  postulate  about  parallel  lines.  See 
text. 

A  recurrent  orbit  in  L  has  length  c  —  irl(er/l  —  e~r/l)  in  which  r  controls  the  intrinsic  metric  of 
the  manifold  as  in  Fig.  4  and  /  represents  an  extrinsic  unit  of  length;  c  grows  in  a  fractionally  exponen¬ 
tial  manner  with  r.  Since  e’1'  —  1  +  r/l  +  1/2 (r/ / ) 2  +...  and  e~r/l  -  1  -  r/l  +  1/2 (r/l)2  -..., 

/  —  2nr(l  +  1/6 (r1  /  I1)  +  ...),  /  -  2 rrr  only  for  a  small  ratio,  r/l.  The  geometric  divergence  of  the 
metric  of  L  "linearizes"  the  measure  on  the  geodesics  on  this  manifold,  the  characteristic  behavior  of 
orbits  of  hyperbolic  sets  (such  as  Smale’s  Anosov-like,  non- wandering  sets),  via  fractional  linear 
transformation,  T .  We  note  that  4>  —  0  and  R 2  — *  L  as  r/  /  — *  °° . 

A  hierarchy  of  recurrent  orbit  lengths,  c  —  f(r)  in  L  space  can  be  visualized  on  the  manifold  of 
Beltrami’s  pseudosphere,  Fig.  5,  the  first  geometric  realization  of  Lobachevskii’s  abstract  fifth  axiom.  It 
was  constructed  using  a  tractrix,  a  curve  whose  diameter  £  —  |x  -  —  re~a,  is  defined  by  the  equidis¬ 

tance  of  all  its  tangents  to  the  xqt  axis  which  Hilbert  proved  could  not  be  extended  without  encounter¬ 
ing  a  fixed-point  singularity,15  in  a  time-reversed  dynamical  contraction  analogous  to  the  expanding 
dynamics  around  the  thermodynamic  critical  point  portrayed  in  Fig.  2. 


Fig.  5  The  normalization  of  length,  f ,  in  Lobachevskii  space  as  realized  on  Beltrami's  pseudosphere.  See  text. 
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The  intrinsic  metric  of  this  geometric  object  representing  L  space  can  be  demonstrated  using  a 
fractional  linear  transformation  T  which  moves  across  diameter  £  carrying  circles  into  circles  and  lines 
into  lines  in  a  smooth  series  of  geodesic  lengths  going  as  +  a)/(l  +  a£„),\a\  <  1  and 

dependent  upon  K  <  0  of  the  hyperbolic  plane.  The  r-dependent  series,  f0-  a; \\  —  2a/l+a2; 
f  2  -  ft  +  a/(l  +  a);f  —  (3a  +  a3)/  (1  +  3a2)...  generates  points  that  converge  to  the  endpoint  of 
the  diameter  as  they  recede  to  infinity  (Fig.  5b) 

Escher’s  Circle  Limit  /  (Fig.  6  left)  portrays  the  geometrically  proportional  rescaling  of  the  hyper¬ 
bolic  plane  in  L  space  when  viewed  on  R 2  as  the  ratio  r/l  is  extended  from  its  Euclidean  neighborhood 
at  r0.  L-1  space  is  seen  in  Escher’s  Smaller  and  Smaller  (Fig.  6  right)  and  portrays  the  opposite  kind  of 
scaling  series.  Again  we  note  that  these  isometries  of  L  would  generate  the  appearance  of  physically 
promising  metric  universality  if  treated  with  renormalization  group  transformations  on  R2. 


Fig.  6  Escher’s  Circle  Limit  (left)  and  Smaller  and  Smaller  (right)  portraying  the  intrinsically  renormalizing 
actions  of  L  space. 

In  hyperbolic  L  space  as  a  one-dimensional  complex  manifold  with  non-Euclidean  metric,  charac¬ 
teristic  of  isotopic  systems  at  their  critical  points  (Fig.  1),  conformal  invariance  (angles  between  any 
two  curves  not  at  vertices  remaining  unchanged)  is  demonstrated  using  fractional  linear  transformation 
of  the  form  T (z)  =*  az  +  b/bz  4-  a  with  a, a  and  b,b  representing  complex  conjugate  pairs,  and 
M2  -  |£|2  =“  1.  The  relation  is  more  easily  visualized  as  T  (r)  =»  e'°  ( z  -  z0)/  (1  —  7(f)  showing  dilation 
as  well  as  rotation  a,  the  angular  derivative  arg(z').  The  distance,  d,  on  a  disc  in  L  space  from  0  to  r 
>  0  using  a  1:1  conformal  mapping  function,  f:L—  L,  and  Caratheodory’s  classic  result  that  |Fr(0)/F(0)j 
<  2  for  all  geodesic  radii  (Fig.  7)  graphed  here  as  circles  (horocycles)  orthogonal  to  (zj  —  1,  is 

r 

J {2dr/ (1  -  r2))  or,  more  familiarly,  log(l  +  r/l  -  r).16  This  is  analogous  to  the  Poincare  metric 
o 

which  reduces  distances  s(r)  on  L  space  like  ds  =  |d4/l  -  |zj2  on  a  manifold  of  K  —  —1.  The  renormal¬ 
ization  of  geodesics  on  this  space  of  constant  negative  curvature  (—  Anosov  flows)  is  intrinsic  to  the 
metric. 


Fig.  7  Horocycle  cascade  on  the  Lobaschevskii  plane  composed  of  concentrically  scaling  circular  geodesic  arcs. 
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Differentiable  Embedding,  Iterative  Roots,  and 
Mode-Scaling  Universality  in  Allosteric  Proteins 

Universality  in  L  space  can  be  sought  in  the  form  of  conformal  invariants  using  their  extrema! 
properties.  If  f.L  — *  L  defines  a  homeomorphic  conformal  map  of  the  unit  disc  to  the  open  disc  which 
is  normalized  as  in  the  above,  !f(0)j  -  1,  there  is  a  least  upper  bound  (and  greatest  lower  bound)  of  the 
radii,  r,  ,ri+1  —  ne~a,  which  has  been  estimated  by  Bloch17  to  be  in  the  range  0.5666  -  0.65647  in  order 
for  distinct  transformation,  f, ,  to  achieve  distinct  values  at  their  centers  (univalence).  Viewing  f:L— *  L 
as  a  shift  transformation  in  a  one-dimensional  dynamical  system  on  the  complex  hyperbolic  plane,  con¬ 
tact  can  be  made  with  comparable  quantitative  universality  at  analogous  thresholds  for  the  change  in 
character  of  topological  conjugacy  functions  at  this  boundary  in  other  one-dimensional  dynamical  sys¬ 
tems. 

Our  theory  is  that  the  phase  transition  of  protein  dynamics  from  resting  sensitivity  to  coherent 
action  consists  of  a  reversible  embedding  of  a  hierarchically  discrete  into  a  continuous,  one-parameter 
group  of  transformations,  C°— ■ -  C=* 1  in  R  and  C,  f*()  F,  (•,  n).  We  use  the  criteria  of  the 
existence  of  the  mth  iterative  root  g  (for  every  natural  number  m,  m  >  1)  for  the  mapping  of  g  such 
that  gm  —  f  for  this  embeddability  in  order  to  reach  another  estimate  of  mode  scaling  exponent  a.18 
Except  for  the  identity  gm  —  /,  the  existence  of  iterative  roots  is  rare,  so  the  argument  requires  the 
assumption  of  the  evolutionarily  selective  restriction  to  a  unique  neighborhood  in  parameter  space  more 
characteristic  of  biological  dynamics  and  not  the  genericity  implicit  in  a  theorem.  For  example,  the 
iterates  of  g  are  the  functions  gn  such  that  gn+ 1  —  gQg”  and  for  m  <  2,  g  is  an  iterative  root  of  order 
m  of  g  if  gm  “  f.  Using  combinatorial  methods,  it  has  been  proved  that  the  quadratic  polynomial 
defined  on  the  complex  plane,  g(g(z))  —  uz2  +  bz  +  c,  has  no  iterative  roots19  and  therefore  cannot  be 
the  form  of  the  resting  state  with  the  discrete-continuous  embedding  property  required  by  allosteric 
proteins. 

In  a  parameter  range  of  measure  zero,  a  contrasting  quadratic  transformation  on  the  real  line 
demonstrates  the  iterative  roots  (factorization)  required  for  embeddability:  geometrically  proportional 
fractional  scaling  of  the  expanding  modes,  the  convolutionarily  self-similar  series  —  <t>',+  1...  such 
that  cj)"-1  +  <t>"  •»  Iterative  roots,  i.e.  (acyclic)  solutions  to  Abel’s  functional  equation 

</>[fOr)l  *  <//(*)  +  l,1*  exist  for  x  -  <J>  because  i/>(<J>2)  -  i M<t>)  +  1,  4>3  -  «t>2  +  <J>  ... 

Beyond  values  for  r  that  generate  subharmonic  bifurcations,  the  quadratic  map  on  the  real  unit 
interval,  f,:R—  R,  [0,ll,  f,(x)  —  rx(l  —  jc),  reveals  its  first  periodic  points  not  a  power  of  two  at  r  ^ 
3. 57. 20  In  parameter  regions  3.56  <  r  <  4  there  is  an  infinite  number  of  periodic  points  along  with 
loss  of  diffeomorphic  conjugacy  (Sarkovskii  Theorem,  1964).  At  r  —  3.83  numerical  analyses  reveal  a 
period  three  whose  points  are  distributed  as  an  isosceles  triangle  (72°,  72°,  36°)  on  the  quadratic  mani¬ 
fold  which,  after  complexification  on  the  half-plane  lm  Z  >  0  and  fractional  linear  transformation  T 
onto  the  boundary  of  the  unit  circle,  the  disc  in  the  Poincar6  metric,  is  an  isometry  equivalent  to  a  rigid 
rotation  (Fig.  8a).  The  sides  are  in  ratio  AC/AB  —  (V 5- 1 )/ 2  —  0.618... 


A  =  oo 


Fig.  8  Period  three  of  the  quadratic  map  on  the  unit  interval  distributes  like  an  isosceles  triangle  on  the  unit 
circle  with  the  Fibonacci  relation  between  longest  and  shortest  sides  (a);  periodic  points  as  the  trace  of  the 
transition  matrix  collect  as  a  Fibonacci  series  (b).  See  text. 
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As  a  subshift  of  the  shift  map  on  two  symbols  at  r  —  3.83  (Fig.  8b),  symbolic  dynamics  show 
that  regions  b  and  c  represent  a  non-wandering,  closed  invariant  set.  When  A  *  |?  |j,  representing 

the  b  and  c  transitions,  is  exponentiated,  we  have  |  2)’ (2  3)'  (3  5) >  (5  g] ^  (g  13)  —  Tr(A),  the 
sequence  of  sums  of  the  diagonal  entries,  collects  finite  forward  orbits  that  grow  in  length  in  the 
Fibonacci  recursive  relation  of  geometrically  proportional  numbers  converging  to  0.618...  *  4>-1.21  The 
transformations  (g  -  <t>"  — 1 ►  <t',+  I)m  — *f  are  self  mappings  g  such  that  gm  —  f  and  g  is  the  mth  iterative 

root  which  exists  for  every  natural  number  m  as  in  a  sequence  of  geometric  means,  M 


KH 


This  factorability  defines  the  capacity  for  embedding,  C°— 1 ►  C*!,  discussed  in  the  context  of  the  dynam¬ 
ics  of  allosteric  proteins 


Consonant  with  the  location  of  the  fixed  points  on  the  boundary  (Fig.  8a),  the  invariant  open  set 
complementary  to  the  closure  of  the  forward  orbit  of  these  critical  points  is  an  expanding  dynamical  sys¬ 
tem.22  Remove  these  three  fixed  points  from  the  invariant  manifold,  and  everything  between  them 
becomes  negatively  curved  into  the  natural  hyperbolic  Lobachevskii  space  of  the  C^1  functioning  sys¬ 
tem  (Fig.  8a).23  The  embeddable  hierarchy  with  iterative  roots  is  C°,  the  discrete  signal-sensitive  rest¬ 
ing  condition,  represented  by  Fibonacci  rotations  <t>"  on  the  boundary,  and  the  expanding  hyperbolic, 
structurally  stable  system  within  the  domain  of  the  non-Euclidean  metric  is  the  manifold  of  transiently 
coherent  protein  action. 


A  Generalization  of  the  Unique  Geometric  Universality 
of  the  C°,  C> 1  Boundary 

(a)  From  the  point  of  view  of  ergodic  theory,  the  unit  circle  is  the  same  as  the  unit  interval.  f:R 
—  R  defined  by  fux  «■*  x  +  w  (mod  1)  is  the  fractional  part  of  x  +  a>  on  10,1]  as  it  is  in  a  rotation  of 
angle#  on  z  €  C,[z|  =*  1,  f:S'  —  S'  defined  by  ( J9  =  e‘2ir(e+u}.  If&>  is  rational  fM  is  finitely  periodic;  if 
c o  is  irrational  the  point  set  resulting  from  the  iterative  composition  of  expanding  map  f  covers  the 
whole  line.  The  two-parameter  nonlinear  circle  map  with  a  quadratic  maximum,  f^S1  —  S1,  defined 
by  f,M  ~  G  +  a  —  r/2n  sin  2n9,  locates  the  boundary  of  the  isometries  defined  by  fractional  linear 
transformations  T  of  L  space.  Poincare  defined  the  rotation  number  for  a  homeomorphic  map  on  the 
circle  S1  asp(frju)  —  lim  (f"  (0)  -  9)/n,  and  at  the  limit  of  topological  conjugacy,  r  —  1  andp(f)  «= 

0.61 8... 24  25  26  The  power  spectrum  of  the  critical  map,  the  "universal  spectrum"  at  r  =  1,  scales  like 
- - a_IS(aaj),a  “  0.618... 


The  unit  circle  can  be  suspended  as  a  torus,  and  then  dividing  the  torus  by  an  involution  f (r )  — * 
— Z  (T2—  S2—  D2)  yields  a  hyperbolic  L-manifold  with  metric  scaling  properties  as  a  function  of 
|rc  —  r|  (Figs.  3  -  6).27  At  high  nonlinearity  the  crinkled  circle  of  numerical  studies  of  critical  maps  on 
R2  demonstrates  similar  hierarchies  before  loss  of  topological  conjugacy  2g.  A  quantitative  similarity 
between  the  Bloch  constant  estimate  on  the  limit  of  univalence  of  one-dimensional  dynamics  on  L  and 
topological  conjugacy  of  the  nonlinear  circle  map  on  the  boundary  of  L  space  is  suggested.  The  Jordan 
curve  resulting  from  the  limits  of  the  e -perturbation  z  —  z2  +  e  (Fig.  9)  is  the  same  quasicircle 
observed  with  critical  nonlinearity  of  fr<u.  With  p  (f)  irrational,  the  dynamics  compose  an  invanant 
minimal  set,  perfect  and  homogeneous,  which  is  either  S1  or  a  Cantor  set  at  r  ^  l.29  The  irrational 
p(D  of  "the  last  torus  to  go"  inscribing  the  universal  spectrum  with  scaling  exponent  (winding  number) 
a  —  0.618...  and  continued  fraction  expansion  [1,1, 1,1...]  can  be  related  to  the  Bloch  constant  limit  on 
the  modular  radii  of  univalent  transformations  T  via  the  well  known  relation  between  T(z)  — 
az  +  b/cz  +  d  and  continued  fractions  [...]  which  can  be  used  to  expand  p(f):  if p(f),  — *  [n , ,  n2  •  •  ■  ] 
andp(0;  — *  [ot|,  m2  •  •  •  ]  thenp(f)^  —  rp(f),  if  there  exist  k,l  such  that  (— 1)*+'  —  1  and  nk+,  and 
m/+,,  for  r  >  0.30 


Fig.  9  When  rJ  —  z1  +  t  the  unit  circle  becomes  an  unrectifiable  Jordan  curve  whose  Hausdorfif  dimension 
varies  continuously  as  « 2  53 


-8- 


(b)  The  one-dimensional  dynamical  system  on  complex  space  comparable  to  the  logistic  map  on 
the  reals  f,:C  — ►  C,  f,  (z)  ■■  rr  (1  —  2),  or  its  conformally  equivalent  z  — 1 •  z2  +  e,  manifests  the  same 
renormal:zable  geometry  on  R2.31  Its  expanding  property  with  respect  to  both  hyperbolic  stretching  of 
the  manifold  (Fig.  3b)  and  the  role  of  critical  loss  of  rectifiability  can  be  seen  in  comparing  the 
undecomposable  continuum  of  orbits  of  the  unstable  homeomorphism  f(z)  —  z2,32  with  absolute 
derivative  two  and  topological  dimension  one  with  its  perturbation  to  a  Jordan  curve  (Fig.  9a, b)  by 
f(z)  —  z2  +  «.  As  e  grows  further,  the  set  has  topological  dimen$:on  one,  linear  Lebesgue  measure 
zero,  and  positive  Hausdorff  dimension  dH,  indicating  that  the  density  scales  as  radius  r°  with  dH  —  1 
+  «2.33  The  Hausdorff  measure  transforms  analogously  to  the  conformally  equivalent  Lebesgue  meas¬ 
ure,  fi,  f*/n  —  |  f|2M  in  the  form  f V  —  |  f(*/Li-  Recent  numerical  simulations  of  the  interior  of  the 
Mandelbrot  bifurcation  set  of  z  — *  z2  +  e  (H-O.  Peitgen,  personal  communication,  1986)  demonstrate 
that  in  going  to  the  boundary  dH  ranges  from  1  to  2  (in  parameter  space)  and  loses  topological  conju- 
gacy  at  a  radius  with  respect  to  the  boundary  which,  estimated  from  the  computer  graphic  display,  is  — 
0.62.  This  may  represent  a  numerical  way  to  sharpen  the  range  of  the  Bloch  constant  estimate  on 
univalence. 

(c)  The  earliest  demonstration  of  the  geometric  origin  of  the  universal  boundary  between  C° 
hierarchical  complexity  and  analytic  expansions  (or  contractions)  on  C  and  therefore  on  R2  was  the 
first  proof  of  the  boundary  of  disjointness  as  a  mode-scaling  invariant  in  C  orR2  (0.618. ..2/3)  by  Cay¬ 
ley.34  He  used  only  the  constancy  of  the  ratio  of  distances  of  imaginary  points  a,  on  any  circle  to  any 
two  antipoints  n ,n'  on  a  common  indefinite  line  and  the  doubling  of  the  angle  between  them  associated 
with  convergence  in  the  geometric  series,  axn/kaxn',  a 2n/ k2a2n',  a3n/k4a3n'...  to  prove  a  sufficient  but 
not  necessary  condition  on  the  "bad  set"  in  the  modulus  of  the  difference  between  two  complex  quanti¬ 
ties,  mod.  (z,  —  Zj)  <  2/3  mod.  zy  ,  for  the  regular  convergence  "of  any  initial  point  taken  at  pleasure," 

2,  —  Zj ,  solving  (z,  +  h)2  —  z2  as  h  — >  0  by  iterating  the  Newton-Fourier  algorithm  z{  — 

(z,2  +  zy2)/2z, ,  z{...  He  showed  that  with  initial  values  in  the  "unfit  segment,"  long  irregularly  behaving 
transients  could  (but  would  not  necessarily)  occur  before  the  eventual  convergence  to  the  attractor’s 
fixed  point.  This  appears  to  be  one  of  the  first  descriptions  in  the  literature  of  the  fractal  set  at  the 
mutual  boundary  of  basins  of  attraction.35-36  In  this  context  the  quantitative  issue  becomes  one  of 
estimating  the  normalized  size  of  the  conformally  invariant,  non-complicated  domain  between  the 
basins  of  attraction  of  the  roots  of  a  complex  quadratic  polynomial,  the  normal  set  of  Fatou.  The  impli¬ 
cation  is  that  the  unstable  fixed  points  (—  the  Julia  sets)  which  generate  fractal  basin  boundaries  analo¬ 
gous  to  those  found,  for  example,  in  nonlinear  second-order  maps  on  the  reals,37  are  complementary  to 
expansive  transformation  groups  with  emergent  geometric  rigidity  such  as  invariant  dH  occurring  at  a 
critical  distance  from  the  Euclidean  approximation  as  the  forward  orbits  escape  to  infinity  (Figs.  4  and 
5). 38  The  unstable  critical  points  of  the  Fibonacci  recursion  seen  in  the  logistic  map  at  r  —  3.83  in  R 
(Fig.  8b)  can  be  identified  with  the  closure  of  the  repelling  periodic  points  of  the  Julia  set  in  C. 

(d)  Another  way  of  generating  hierarchies  at  the  C°,  C^1  boundary  includes  the  Hardy- 
Weierstrass  trigonometric  map,  f:S'  —  S1,  defined  by  -  I"  cos  (A”r)  with  the  Holder  exponent 
serving  as  the  bound  on  differentiability,  logy-1/log  A  »»  l.39  We  note  here  that  (0.618. ..)“'  = 

1.618...  uniquely.  A  generalization  of  this  lacunary  power  series,  a  transformation  of  the  cycles  of  all 

00 

orders  analogous  to  that  seen  in  the  logistics  map  at  r  *  3. 83, 40  21a, z"'  i/n,  >  k  >  1)  with 

1 

equivalence  of  its  Hausdorff  and  spectral  dimensionalities  has  been  called  a  perfect  set  of  unicity  and 
multiplicity.41  We  shall  see  below  that  this  relation  is  realized  in  physical  studies  of  proteins.  A  critical 
level  of  competition  between  convergent,  y,  and  divergent.  A,  motions  in  these  trigonometric  maps  of 
hyperbolic  orbits  in  thermodynamic  phase  space  represents  the  forces  of  hydrogen  bond  maintenance 
versus  solvent  entropy  preservation  and  polymeric  aggregation  versus  chain  entropy  maximization  con- 
s’itming  the  hyperbolic  pressure  of  allosteric  proteins  in  solution.  This  dissipative  yet  area-preserving 
dynamic  in  R2  can  be  represented  most  simply  by  a  hyperbolic  linear  map  fOr)  -  Ax,  A  -  |q  ,  0  < 
y  <  1  <  A,  and  det  A  —  1.  A  single  parameter,  derived  from  the  sequence  of  amino  acids  in  peptides 
as  hydrophobicities,  determines  the  hyperbolic  pressure  of  proteins  defined  asr  -  <]  A  |  —  |  y  |)  >  0. 

(A  +  y )  =»  0  (see  below) . 


A  simple  construction  (Fig.  10)  generalizes  the  role  of  the  critical  balance  of  the  components  of 
hyperbolic  pressure  around  X  +  y  —  0  as  the  expanding  and  contracting  manifolds  of  a  one- parameter 
group  of  transformations  in  R2  orC,  which  is  analogous  to  the  steepness  (or  height)  of  the  quadratic 
manifold  of  Fig.  8b.  It  demonstrates  the  geometrically  invariant  self-similarity  of  the  critical  boundary 
between  the  C°  and  C*1  global  (not  local)  conditions.  In  the  near-critical  case  when  X  +  y  <  0,  Fig.  10 
shows  three  lines  with  a  diminishing  series  of  segments  with  ratio  1  /<t>  that  have  the  unit  of  measure  as 
the  first  term:  AB  -  BC  +  CD,  BC  -  CD  +  DE...and  AB/BC  -  BC/CD  -  CD/DE...  -  <t>.  fV*2 
—  R2,  leads  to  a  series  of  lines  at  equal  intervals,  0,  (like  time)  radiating  from  the  pole  (fixed  point). 


Fixed 

o 


Fig.  10  A  one-parameter  group  of  transformations,  generates  Fibonacci  scaling  at  the  boundary  of 
differentiability.  See  text. 


The  self-similarities  intrinsic  to  geometric  progression  abound.  The  triangles  PAA'  and  PA'A"  are 
similar,  since  the  corresponding  angles  are  equal.  On  these  grounds,  PA'/PA  =  PA'VPA'...,  and  the 
lengths  PA,  PA',  PA”  constitute  a  geometric  series.  The  triangles  PBB',  PB'B"...  are  part  of  the  same 
series  of  triangles,  and  PB,  PB\  PB"  compose  terms  of  the  same  geometrical  progression.  Of  course 
PA,  PB,  PC...  form  a  geometrical  progression,  for  they  are  "equally  spaced*  terms  of  the  one  just  men¬ 
tioned  (see  Fig.  4).  Another  geometrical  progression  is  formed  by  the  lengths  AA',  A'A",...BB\  B'B”... 
The  quadrilaterals  AA'B'B,  A' A"  B"B'...  are  all  similar  since  they  have  equal  corresponding  angles  and 
pairs  of  sides  in  proportion. 

Although  this  construction  begins  with  the  one  parameter  being  the  perpendicularity  of  the  A—* 
A',  A'—  A"  lines  intersecting  with  the  radii,  a  Fibonacci  recursion  of  radial  lengths  would  have 
emerged  naturally  from  P  if  the  parameter  (instead  of  sequential  90°  interactions)  was  the  restraint  that 
PA,  PA',  PA”  be  in  geometric  progression.  In  the  same  way,  the  problem  of  determining  an  infinite 
number  of  countable  intermediate  points  reduces  to  that  of  finding  the  geometric  means,  M  =* 

tt  i/m 

Y[ak  ,  between  the  lengths  of  existing  radii  PA, PA'...  These  one-parameter  homothetic  transfor- 

k 

mations  put  the  iterative-root,  embeddability  theorem18  relating  the  geometrically  singular  condition  on 
’he  common  boundary  between  discreteness  and  continuity  in  expanding  (or  contracting)  maps  of  the 
plane  in  the  context  of  Ahlfors’s  "ultrahyperbolic”  Poincare  metric  in  L  space  with  Gaussian  curvature 
K  —  -1  via  conformal  invariance.27  The  elements  of  the  group  T(z)  C  L  are  conjugate  to 
f*9o«(x  y)  C  R2  on  the  geometric  space  of  Fig.  10.  Both  contain  the  action  of  renormalization  as  an 
intrinsic  feature  of  the  geometry  of  their  manifolds;  both  are  transformationally  invariant  with  respect 
to  translation,  inversion,  rotation,  and  dilation. 

If  in  Fig.  10  the  sequence  of  radii  from  circle  center  P,  PA,  PA',  PA"...  were  not  in  geometric 
ratio,  but  PA/PA'  —  PA'/PA"  —  1,  the  triangles  PAA'  and  PA'A”  as  well  as  their  nested,  self-similar 
x-uiivergence  PBB',  PB'B",  PCC,  PC'C"...  would  satisfy  the  strong  triangle  inequality  d(PA)  —  max  d 
(PA')  ^  AA’).  This  generates  a  sequence  of  isosceles  triangles  that  serve  as  a  condition  on  the  boun¬ 
dary  of  individual  clusters  in  hierarchical  structures  as  defined  by  relative  nearness  on  a  topological 
space.42  This  boundary,  then,  defines  the  limit  of  disjointness  in  ultrametric  spaces,  hierarchical  lattices 
like  corporate  tables  of  organization  or  upside-down  trees.  This  topology  with  the  common  currency  of 
time  has  been  used  to  relate  the  stretched  exponential  relaxation  of  hierarchical  dynamics  (the  charac¬ 
teristic  function  of  the  fractional  power  law  scaling  spectra  of  L6vy,  convolutionarily  self-similar  distri- 
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bution  laws43)  and  random  walks  on  hierarchical  manifolds.44  43  Isomorphisms  of  the  L6vy  mode  distri¬ 
bution  are  generated  by  the  Hardy- Weierstrass  cosine  series  39  as  well  as  the  Salem  sets  of  uniqueness 
and  multiplicity41;  at  their  boundaries,  both  demonstrate  the  nested  property  of  the  strong  triangle  ine¬ 
quality  of  ultrametric  spaces  and  the  iterative  roots  (like  self-similar  Livy  convolutions)  of  embeddable 
maps  at  the  C°,  C* 1  boundary. 

In  passing  we  note  that  in  the  critical  neighborhood  of  X  +  y  >  0,  PA/PA'  -  PA' A"  >  1  such 
that  angles  similar  to  B'AA'  —  90°  (rather  than  those  similar  to  B'A'A  —  90°  in  Fig.  10),  a  similar 
geometric  proportionality  with  embeddable  scaling  <J>  bounds  the  radii  from  the  other  side,  i.e.  with 
cords  similar  to  AA'  normalized  to  l,the  radii  are  equal  tod*  and  will  grow  both  additively,  <J>"+1  — 

4>  +  <J>",  and  multiplicatively,  <|)',+  1  —  jn  a  manner  characteristic  of  this  unique  geometric  pro¬ 
gression.  With  the  strong  triangle  inequality  of  ultrametric  spaces,  topological  hierarchies  are  bounded 
on  both  sides  by  the  geometric  scaling  of  the  Fibonacci  recursion  relation.  Since  this  boundary  has 
measure  zero,  the  apparent  diversity  in  descriptors  may  be  regarded  as  identities.  A  related  treatment 
of  the  boundary  of  probabilistically  disjoint  spheres  bounded  by  imaginary  quadratic  irrationality,  which 
can  be  represented  by  periodic  continued  fraction  expansions  and  contains  logarithmic  density  distribu¬ 
tions  of  geodesic  orbits,  exploits  the  hyperbolic  metric  of  Figs.  3  -  6.  Sullivan’s  theorem46  that  the 
universality  of  this  boundary  on  disjointness  comes  from  a  one-dimensional  reduction  fits  both  Hamil¬ 
tonian  dynamics  when  studied  on  a  circle  and  the  three-parameter  Cartwright-Littlewood  equation  when 
examined  on  the  plane  (see  below).  Fig.  10  portrays  the  underlying  geometric  universality  at  the 
border  of  differentiability  in  R2  and  C  that  has  emerged  from  a  variety  of  formalisms,  including  the 
renormalization  group  transformation  of  distances  as  the  critical  point  is  approached.  Whether  this 
measure  zero  boundary  looks  fractal  or  smooth  depends  on  the  point  of  view  of  the  observer. 

A  Differential  Equation  for  the  Thermodynamics 
of  Signal-Sensitive  Proteins 

Through  nonlinear  positive  and  negative  dampening,  the  r -parameter  of  the  Cartwright- 
Littlewood,  van  der  Pol  equation47  regulates  the  global  hyperbolic  mechanism  of  the  flow  at  the  C°,  1 

boundary;  the  critical  value  of  the  normalized  amplitude  parameter  is  b  ^  (d>-1  -  3).  This  global 

property  represents  the  competing  components  of  the  hyperbolic  stability  of  allosteric  proteins  in  water. 
The  system  at  r  >  rc  can  be  regarded  as  an  expanding  dynamical  system,  analogous  to  a  polynomial 
vector  field  in  R2  in  x  and  y  when  studied  at  high  nonlinearity.28  Differentiating  f(x)  = 
x  -  r(l  -  x2)x  +  x  —  rb  cos  Xr  and  rescaling,  y/r  — ■  y,  leads  to  the  pair  of  ordinary  differential  equa¬ 
tions:  x  -  y  —  r  (x  -  x3/3)  and  y  —  —x/r  +  b  cos  \t  whose  behavior  can  be  examined  in  phase  space 
as  well  as  on  the  two-parameter  plane.  The  dimension  x  portrays  hydrophobic  packing,  AH,  in  the 
interval  [0,1]  andy,  the  entropy,  T°AS,  normalized  in  the  same  way.  Fast  hydrophobic  aggregation  fol¬ 
lowed  by  the  slower  rearrangement  of  hydrogen  bond  cages  makes  the  relaxation  oscillations  of  this 
equation  set  realistic.  Because  both  properties  are  tightly  linked  (are,  in  effect,  two  aspects  of  the  same 
thermodynamic  motions),  the  behavior  can  be  viewed  as  x  versus  x  trajectories  in  one-dimensional 
dynamics  on  a  two-dimensional  phase  plane.  The  parameter  X  signifies  the  complex  hydrophobic  mode 
of  a  polypeptide  transmitter  or  ligand,  the  normalized  external  "driving"  signal  is  cos  Xr;  b  serves  as  the 
amplitude  parameter  of  the  driving  force  representing  the  amount  of  transmitter/ligand  message  which 
we  use  in  a  pharmacological  "dose-response"  fashion  to  investigate  the  geometric  universal  scaling  rela¬ 
tion  at  the  boundary  of  differentiable  mode-locking  behavior.  The  physical  exchange  between  polypep¬ 
tide  transmitter  and  protein  receptor  is  in  the  form  of  interfacial  surface  tension  generated  by  their 
mutual  hydrophobic  attractive  forces  which  are  lOOx  greater  than  predicted  by  van  der  Waal  theory  and 
extend  up  to  10  nanometers,  the  size  of  a  neuronal  gap  junction.48 

At  b  —  0,  numerical  analyses  of  this  system  showed  that  increasing  hydrophobic,  hyperbolic  pres¬ 
sure  r  stretches  the  manifold  of  the  unperturbed  relaxation  oscillator  (Fig.  I  la)  and  slows  its  native  fre¬ 
quency  in  a  nonlinear  manner  (Fig.  lib). 

Power  spectral  transformation  as  a  density-of-states  spectrum  over  increasing  hyperbolic  pressure 
parameter  r  is  qualitatively  1/f,  and  a  deconvolution  generates  the  characteristic  stretched  or  algebraic 
correlation  function  of  hierarchical  dynamical  structures.  It  is  of  more  than  historical  interest  that  van 
der  Pol’s  original  diode-rectifier  circuit49  manifested  an  r-dependence  of  the  relaxation  times  r,  that 
varied  as  r,+  1/r,  —  <J>z/<y<f  Here  r/u$  —  RC,  and  emf- driven  RC  circuits  model  the  dynamics  of  ion- 
conductance  behavior  of  allosteric  membrane  proteins  very  well.50 


Fig.  11  An  increase  in  hyperbolic  pressure  parameter  r  stretches  the  manifold  (a)  and  slows  the  frequency  (b) 
in  the  unforced  Cartwright-Littlewood  attractor. 

The  characteristic  influence  of  an  increase  in  hyperbolic  pressure  0.5  <  r  <  1.5  on  the  periodi¬ 
cally  transmitter-  or  ligand-driven  system,  b  >  0,  with  K  held  between  the  Arnol’d  tongues51  of  mode- 
locked  winding  numbers  can  be  seen  in  Fig.  12,  left  versus  right.  As  might  be  anticipated  from  Fig.  11, 
an  increase  of  hyperbolic  pressure  (up  to  some  quadratic  limit)  increases  the  complexity  between 
mode-locked  tongues  in  the  phase  space  of  the  orbital  dynamics  from  which  one  would  predict  an 
increase  in  some  measure  of  complexity  such  as  the  Hausdorff  dimension  dH.  This  leads  to  the  physi¬ 
cally  counterintuitive  prediction  that  resting  (not  mode-locked)  proteins  with  higher  average  amino  acid 
hydrophobicities  will  manifest  higher  Hausdorff  dimensionality  in  time  and  space  (see  below). 


Fig.  12  The  Canwright-Liulewood  attractor  between  mode-locking  parameter  zones  at  lower  (left)  and  higher 
(right)  levels  of  hyperbolic  pressure.  See  text. 

The  strange  dynamics  of  the  periodically  driven  attractor  (Fig.  12b)  are  controlled  by  the  r- 
dependent  linkage  between  x  and  x  in  which  the  velocity  adjusts  itself  to  make  the  bivalent  dampening 
balance  the  forcing  term,  i.e.  area-preserving,  dissipative  forcing.  The  global  hyperbolic  pressure  gen¬ 
erated  by  x0,  x0  repulsion  and  x  >  ±  1,  x  >  2/3  attraction  is  focused  at  the  saddle-sinks  where  x  -  0 
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and  x  *■  ±  1  (Figs.  11  and  12).  They  play  a  determining  role  like  critical  points  in  the  region  of  con¬ 
tinuous  phase  transitions  (Fig.  2)  and  points  with  characteristic  derivative  F(z)  “  0  when  in  the  limit 
sets  of  analytic  maps  such  as  z  — »  z2  +  e  or  the  conformally  equivalent  z  —  ez  (1  —  z)  (Cvitanovic  and 
Myrheim,  unpublished,  1984). 31 

The  trajectories  speed  around  the  x  axes,  slow  critically  at  the  saddle  sinks,  undergo  a  periodic 
"tug-of-war''  between  leaving  and  staying  which  generates  whorls  of  splinter  orbits52  that  grow  in 
number  tangent  to  the  insets,  past  the  saddle  sinks  until  they  get  too  close  to  the  central  repeiler  jr0,  x0 
and  are  driven  to  accelerate  toward  the  next  x  axis  and  begin  again  (Fig.  12b).  That  the  complexity 
arises  from  the  invisible  central  repeiler  and  the  mixmaster-like  actions  of  the  two  saddle-point  repellers 
between  linear  segments  of  the  flow  is  analogous  to  the  dynamics  in  the  neighborhood  of  the  Julia  sets 
in  polynomial  maps  of  the  complex  plane  (and  unstable  fixed  points  at  the  fractal  basin  boundary).36 
Whorls  of  intermittent  "bursts  of  persistence"  have  been  reported  as  generic  in  numerical  studies  of 
critical  nonlinear  maps  of  the  plane.53  This  behavior  can  be  quantified  using  local  escape  rates  from  ini¬ 
tial  conditions  in  the  neighborhoods  of  these  fixed  point  repellers.54  Physically,  these  critical  points 
represent  the  extremes  of  the  protein’s  breathing  motions  between  maximal  hydrophobic  packing  and 
nearly  denaturing  entropy  gain  through  unfolding  of  the  peptide  polymer.3 

Squeezed  between  accelerating  repulsion  at  x0,x0  and  dampening  contraction  at  x  >  ±  1,  x  > 

2/3,  the  dynamics  create  an  annulus  on  the  plane,  which  although  dissipative,  is  analogous  to  the  area¬ 
preserving  Birkhoff  attractor  with  different  inner  and  outer  rotation  numbers,  non-trivial  intervals 
between  them,  and  two  distinct  invariant  points.55  The  existence  of  infinitely  many  periodic  motions 
near  a  stable  periodic  motion  as  well  as  motions  that  are  not  periodic  but  are  uniform  limits  on  periodic 
motions  follows  from  the  PoincarS-Birkhoff  theorem  and  is  seen  directly  in  a  variety  of  these 
Cartwright-Littlewood  attractors  (Fig.  12b).  Recent  numerical  studies  of  a  two-parameter  family  of  dis¬ 
sipative  twist  maps  of  the  annulus  (Aubry-Mather  sets)  demonstrated  non-trivial  intervals  enclosed  by 
the  internal  and  external  rotation  numbers  similar  to  those  seen  in  the  conservative  case.56  Their 
existence  in  these  innately  unstable  analog  computer  simulations  of  the  Cartwright-Littlewood  attractor 
confirms  previous  observations  that  these  rotation  intervals,  as  studied  in  digital  simulations  of  dissipa¬ 
tive  circle  maps,  are  remarkably  stable  to  € -perturbation57  as  well  as  to  random  noise.  This  may  justify 
referring  to  Fig.  12b  as  a  phase-locked  strange  attractor  (M.  Casdagli,  personal  communication,  1986). 

It  portrays  the  complex  thermodynamic  phase-space  dynamics  of  the  hydrophobic  free  energy  of 
signal-sensitive  proteins  where  AG  ^  0. 

The  analogy  between  the  behavior  of  this  f(x)  in  the  intervals  of  rotation  numbers  on  R2  and  that 
ofT(z)  on  L  is  clear.  ForT(z)  -  az  +  b/  {cz  +  d)  with  ad  -  be  =  1  andr(z)’-  (cr  +  d)~2,  the 
orbit  (zjcr  +  d|=  1)  is  the  isometric  circle,  c,  ofT(z)  since  T(z)  rescales  lengths  within  it.27  Or,  from 
another  point  of  view,  as  in  Figs.  4  -  7,  fr  it2  — •  L  is  a  transition  from  a  Euclidean  to  a  Poincare  metric 
at  some  critical  hyperbolic  pressure  r.  This  R2  — *  L  region  of  transition  is  the  same  as  that  of  degen¬ 
erate  Hopf  bifurcations.58  In  the  L  neighborhood,  with  a  surface  of  constant  negative  curvature,  the 
orbits  have  positive  topological  entropy  and  the  system  demonstrates  near  Anosov  structural  stability,59 
I  behavior  on  an  annulus  analogous  to  that  of  a  hyperbolic  automorphism  on  the  torus.60  The  implicit 

1  analogy  with  Axiom  A  flows  as  the  generic  hyperbolic  sets  of  allosteric  proteins  can  be  argued  more 

easily  if  we  weaken  the  conditions  so  that  no  more  than  one  characteristic  exponent  is  equal  to  zero, 
leaving  exponential  growth  and  decay  of  perturbations  intact  anu  eliminating  the  strong  uniformity  con¬ 
dition  which  is  physically  difficult  to  realize.  One  can  visualize  the  center  manifold  of  protein  actions 
such  as  catalysis  or  binding  modulated  by  transverse  hyperbolic  macromolecular  motions.3  In  the 
region  of  geometric  proportionality,  distance  is  rescaled  via  the  Fibonacci  recursion  relation  (Fig.  10). 
This  neighborhood  of  the  C°,  C* 1  boundary  is  that  of  the  signal-sensitive  condition  manifested  in  the 
universal  power  spectrum  of  allosteric  proteins.  With  loss  of  hyperbolic  pressure,  r  <  rc,  fr  :L  —  R2, 

,  the  topological  entropy  h(f)  goes  to  zero  with  the  loss  of  the  negative  curvature  as  L  —  R2.61  A  similar 

[  l°ss  of  complexity  occurs  in  the  (rA<&)  zone  of  mode  locking  on  the  invariant  circle,  rational  rotations 

!  on  the  boundary  c. 

j  Numerical  analyses  of  the  Cartwright-Littlewood  equation  in  r  A  space  demonstrate  at  0  <  r  < 

1.0  almost  the  full  Farey  sequence  of  winding  numbers  p(f),  up  to  a  denominator  of  9,  withp(f),  that 
appear  to  be  continuously  dependent  on  r,  A,  and  b.  Generally,  increasing  r  ^  0.5  decreases  parame¬ 
ter  space  complexity,  seen  in  the  progressively  lower  denominators  of  p(f),  which  increase  in  size  and 
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decrease  in  number  in  a  process  seen  also  in  the  expanding  Amol’d  tongues  representing  the  mode 
locking  zones  of  the  winding  numbers  (1/1,  2/1,  3/1...  going  from  left  to  right  in  Fig.  13.  With  the 
exception  of  the  1:1  tongue,  the  continuous  parameter  dependence  of  p(f)  —  f(r,X,6)  appears  to 
break  down  into  islands  across  X  as  a  function  of  increasing  r  and  b . 


Fig.  13  Mode-locking  zones  of  the  Cartwright-Littlewood  attractor  in  r  x  parameter  space.  Zones  of  1:1.  2:1. 
3:1,  4:1...  are  indicated  in  a  dose-response  study  of  the  system  over  increasing  b  The  1:1  Arnol'd  tongue  goes 
across  all  of  r  at  0.666  <  b  <  0.700. 


Like  the  existence  of  non-trivial  intervals  between  internal  and  external  rotation  numbers 
(Fig.  12),  at  high  nonlinearity  r  the  loss  of  tongues  and  the  emergence  of  mode-locked  islands  sur¬ 
rounded  by  complicated  sets  (here  overlapping  of  the  neighborhoods  of  at  least  twop(f),)  are  reminis¬ 
cent  of  the  behavior  of  Hamiltonian  dynamics  with  perturbation-induced  destruction  of  KAM  surfaces 
and  the  appearance  of  islands  surrounded  by  Amol’d  diffusive  stochasticity  in  the  neighborhood  of  their 
separatrices.25  A  third  similarity  to  the  results  of  studies  of  nonlinear  stability  in  conservative  systems 
is  the  scaling  relation  of  the  invariant  circles  at  the  boundary  of  topological  conjugacy  (previous  sec¬ 
tion).  We  see  a  transition  to  1:1  mode  locking  across  all  of  r  in  the  general  vicinity  of  the  driving  fre¬ 
quency,  amplitude  scaling  parameter  b  —  2/3  as  first  reported  by  Cartwright  and  Littlewood.47  This 
scaling,  comparable  to  the  mode  ratios  and  power  law  (<t>,  <t-1)  generated  by  a  Fibonacci  recursion  at 
the  boundary  of  chaos  in  standard  map  models  of  KAM  systems,25  has  been  noted  previously.  Numeri¬ 
cal  simulation  using  double  precision  digital  computing  will  be  required  to  determine  the  boundary  pre¬ 
cisely.  The  mode-locked  insensitivity  at  b  >  2/3  models  the  habituation-desensitization  behavior  of 
:gr.al-sensitive  systems  at  high  levels  of  transmitters,  ligands,  or  drugs.2 

A  link  to  the  low-dimensionaJ,  thermodynamic  physics  of  proteins  in  aqueous  solution  can  be 
made  from  the  geometric  universality  of  (<t> ...2/3)  on  R2  and  C  seen  in  both  dissipatively  forced  and 
conservative  nonlinear  systems  to  the  enthalpy-entropy  compensation  of  critically  forced  water  by  recal¬ 
ling  that  the  mode  ratio  w/ai]  found  in  the  Fourier  transformed  velocity  measurements  made  in  the 
Couette-Taylor  system  (cii2  was  transient)  at  r/rc  «■  10  -  11,  just  before  the  appearance  of  the  chaotic 
band,  is  in  the  range  of  0.62  -  0.65.62 


Stapleton’s  Hausdorff  Measure  and  Spectral  Scaling  Exponent 
on  Allosteric  Proteins  in  R2  andC 

The  preceding  development  leads  to  four  rather  straightforward  predictions63-64:  (1)  The  tempera¬ 
ture  dependence  of  the  mode  scaling  exponent  of  the  density-of-states  spectrum  of  allosteric  proteins 
will  reflect  their  hyperbolic,  hydrophobic  pressure  at  the  limit  of  the  C°,  C* 1  boundary  as  a  Fibonacci 
recursion  relation,  (2)  the  abstract  geometric  grounds  of  these  relational  constants  will  lead  to  experi¬ 
mentally  verifiable  isometries  between  R2  andC,  (3)  the  physical  basis  of  this  relational  universality  is 
the  hyperbolic  thermodynamic  stability  of  allosteric  proteins  in  aqueous  solution  which  will  be  quantita¬ 
tively  related  to  their  complex  hydrophobicities;  and  (4)  these  measures  on  the  complex  hydrophobic 
dynamics  of  proteins  will  be  reflected  in  (all)  other  proper  indices  of  their  thermodynamic  structure.10 

The  temperature  dependence  of  the  density-of-states  spectra  in  studies  of  Raman  relaxation  rates 
in  C  of  four  iron-binding  allosteric  proteins  when  compared  with  ordinary  solids  was  found  to  be 
anomalous.65  The  density  varied  with  temperature  past  both  the  T2  and  characteristic  phonon  localiza¬ 
tion  regions  of  the  scattering  relation  of  amorphous  solids  asy™-1,  and  m  varied  from  1.61  ±  0.05  to 
1.67  ±  0.03  (i.e.  a  complex  temperature  scaling  exponent  a  of  0.618. ..2/3).  The  Hausdorff  dimension 
dH  of  the  x-ray  crystallographic  maps  of  the  distributions  of  the  amino  acid  alpha  carbons  of  these  same 
proteins  in  R2  determined  via  packing  density  as  a  function  of  the  radius,  mass  —  r“,  yielded  nearly 
identical  values  fora.66 

The  hyperbolic  thermodynamic  stability  hypothesis  predicts  that  in  contrast  to  the  intuition  that 
amino  acid  chains  that  are  hydrophobic  aggregate  to  occupy  smaller  areas  of  Fourier  and  Euclidean 
space,  the  spectral  scaling  and  Hausdorff  mass  scaling  exponent  a  would  vary  with  the  average  hydro- 
phobicity  per  residue  using  a  Tanford-derived  set  of  amino  acid  hydrophobic  free  energies  determined 
by  their  ethanol-water  partitions.67  Table  1  lists  six  iron-binding  proteins  for  which  the  amino  acid 
sequences  were  available.  We  note  that  whereas  allosteric  proteins  such  as  myoglobin  demonstrate 
phase  transitions  and  the<t>  scaling  of  the  C°,  C5*"1  boundary,  ferredoxin  and  putidoredoxin,  which  are 
not  allosteric  proteins,  do  not. 


Protein 

Table  I 

Avg.  Hydrophobicity  (kcal) 
per  Amino  Acid  Residue 

ESR 

(C) 

X-ray 

(R2) 

Myoglobin  azide 

1.34 

1.65 

1.65 

Myoglobin 

-1.34 

1.67 

1.65 

Cytochrome  P-450 

1.32 

1.64 

1.63 

Ferredoxin 

1.11 

1.34 

1.37 

Putidaredoxin 

1.21 

1.34 

1.36 

Cytochrome  C-551 

1.27 

1.43 

1.50 

A  structural  correlate  of  the  complex  hydrophobic,  hyperbolic  pressure  hypothesis  is  that  faster 
frequency  modes  and  increased  spectral  and  Hausdorff  dimension  will  be  reflected  as  well  in  increased 
optical  rotatory  dispersion  as  percent  helix  and  related  in  part  to  the  average  amino  acid  hydrophobicity. 
For  examples:  lysozyme  with  an  average  hydrophobicity  of  1.24  kcal/residue  and  a  dominant  mode 
spectrum  demonstrating  a  2.5-residue  wavelength  (see  below)  has  a  Stapleton  dH  of  1.66  and  35-45% 
helix;  pancreatic  ribonuclease,  with  an  average  hydrophobicity  of  1.02  kcal/residue  and  a  dominant  5.9 
and  11.43-residue  hydrophobic  mode,  has  a  Stapleton  dH  of  1.33  and  less  than  15%  helix.6*  Other 
examples  for  which  data  are  available  on  dH,  %  helix  respectively:  myoglobin  1.65,  87.6%;  adenylate 
kinase  1.59,  61.9%;  cytochrome  C  1.61,  48.5%;  lactic  dehydrogenase  1.50,  41.6%;  carboxypeptidase 
1.51,  39.5%;  subtilisin  1.46,  32.4%;  alcohol  dehydrogenase  1.43,  29.1%;  papain  1.42,  27.8%,  carbonic 
anhydrase  1.33,  16%;  chymotrypsin  1.32,  4.1%;  concana valin  A  1.26,  2.5%. 
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Complex  Hydrophobic  Transformation  of  Amino  Acid 
Sequences  Predicts  Hormone- Receptor  Recognition 

The  dynamics  of  the  Cartwright-Littlewood  equation  for  the  transmitter/ligand-driven,  enthalpy- 
entropy  oscillations  of  allosteric  proteins  with  the  capacity  for  signal-induced  C°,  C* 1  transitions  predict 
that  hydrophobic  mode  resonances  between  polypeptide  hormones  (cos  Xf ,  b)  and  membrane  receptors 
allosterically  regulated  by  r  in  a  common  aqueous  field  are  the  source  of  reversible  action-inducing  and 
(if  too  extreme)  irreversible  desensitizing  phase  transitions. 

Fourier  transformation  of  the  partition  coefficients  as  hydrophobic  free  energies  in  kcal/mol67 
representing  the  amino  acid  sequences  as  a  time  series  (Table2)63,64  generated  hydrophobic  mode  power 
spectra  for  polypeptide  polymers. 


Table  2 

Amino  Acid  Hydrophobicities  (kcal) 


Ala 

0.87 

Glu 

0.67 

Arg 

0.85 

Gin 

0.00 

Asn 

0.09 

Gly 

0.10 

Asp 

0.66 

His 

0.87 

Cys 

1.52 

(le 

3.15 

Leu 

2.17 

Ser 

0.07 

Lys 

1.64 

Thr 

0.07 

Met 

1.67 

Trp 

3.77 

Phe 

2.87 

Tyr 

2.76 

Pro 

2.77 

Val 

1.87 

For  example.  Fig.  14  (left)  shows  the  hydrophobic  sequence  spectrum  of  the  41-amino  acid  cor¬ 
ticotropin  releasing  factor  (CRF)  and  five  analogues  with  —50%  amino  acid  homology  but  nearly 
equivalent  potency69  that  demonstrate  a  common  dominant  mode  of  —  2.1  residues.  Fig.  14  (right) 
shows  similar  transformation  of  the  44-amino  acid  growth  hormone  releasing  factor70  and  several  pep¬ 
tides  with  a  variety  of  different  previously  designated  actions  which,  although  considerably  less  than 
50%  homologous  with  respect  to  specific  amino  acids  are  roughly  equivalent  in  their  capacity  to  release 
growth  hormone  from  pituitary  analog  cell  preparations  (F.  Zeytin,  personal  communication,  1985). 

All  their  spectra  are  dominated  by  4.0-residue  hydrophobic  modes. 
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Fig.  14  Hydrophobic  mode  spectra  of  CRF  (left)  and  GRF  (right)  and  analogues. 


A  hydrophobic  spectral  survey  of  the  neuropeptides  known  to  be  present  in  the  mammalian  ner¬ 
vous  system71  (Table  3)  led  to  the  discovery  of  a  unique  secondary  8.0-residue  mode  in  the  74-residue 
peptide  ubiquitin.  Fig.  15  shows  the  spectrum  of  a  two-point  moving  average  of  the  hydrophobic 
sequence  which  removed  the  omnipresent  —2. 1-residue /3-strand  component  (3.4  -  3.7  represents  the 
a -helical  segments). 


Table  3 

Hydrophobic  Sequence  Wavelengths 

Peptide 

Wavelengths 

Peptide 

Wavelengths 

Primary 

Secondary 

Primary 

Secondary 

ACTH 

2.35 

13.3 

Insulin  A 

4.00 

3.08 

Angiotensin  II 

4.00 

Insulin  B 

5.33 

Atrial  peptin  1 

4.00 

0-LPH 

4.44 

Atrial  peptin  II 

4.44 

LHRH 

2.22 

3.08 

Bombesin 

4.44 

2.22 

Motilin 

2.86 

Bradykinin 

2.22 

3.08 

a-MSH 

2.17 

Calcitonin 

3.64 

0-MSH 

2.67 

CCK-33 

3.08 

Neuropeptide  Y 

3.33 

CRF 

2.11 

4.44 

Pancreatic  polypeptide 

3.33 

2.00 

0-Endorphin 

4.44 

Proctolin 

2.50 

Gastrin 

2.67 

• 

Secretin 

4  44 

3.33 

GH 

3.64 

Somatostatin-28 

13.3 

4.00 

GIP 

5.00 

2.67 

Substance  P 

3  64 

2.67 

Glucagon 

4.00 

13.3 

Ubiquitin 

2.11 

8.00 

GRF 

4.00 

VIP 

4.00 

2.35 

Fig.  15  Hydrophobic  mode  spectrum  of  ubiquitin. 


A  comparable  7  to  9-restdue  mode  was  found  in  the  53-residue  epidermal  growth  factor  (EGF)72 
as  well  as  other  peptides  mitogenic  for  mesodermal  and  ectodermal  tissues  such  as  transforming  growth 
factor  (TGF)  (Fig.  16). 

Only  a  few  receptors  have  been  isolated  and  completely  sequenced;  among  them  is  the  EGF 
receptor  protein.73  The  126-residue  extramembrane  segment  immediately  following  the  24-acid  signal 
sequence  was  found  to  be  dominated  by  a  7  to  9-residue  mode  that  emerged  again  in  the  post¬ 
membrane,  cytoplasmic  domain.  Extensive  homology  of  the  EGF  and  tyrosine  kinase  receptor  related 
to  oncogenesis  suggests  a  potential  role  for  the  7  to  9-residue  mode  in  the  regulation  of  oncogenes.74 

Another  recently  sequenced  receptor,  for  glucocorticoids,75  was  examined  as  a  function  of  the 
average  hydrophobicity  of  each  ribosomal  segment  (—35  residues),  and  strong  7. 5-residue  peaks  were 
found  on  both  sides  of  the  membrane,  along  with  the  common  /3 ,  a ,  and  «  (4.4  -  5)  modes  (Fig.  17) 
This  suggested  that  some  aspects  of  steroid  actions  might  be  mediated  by  a  7  to  9-residue  hydrophobic 
'"ode,  the  implicit  assumption  being  that  the  critical  thermodynamic  behavior  of  water  mediates  the 
mutual  recognition  of  members  of  dissimilar  chemical  families,  here  the  peptide  receptor  and  non¬ 
peptide  hormone. 

We  looked  for  evidence  that  the  7  to  9-residue  mode  that  dominated  ubiquitin,  transforming 
growth  factors  and  their  receptor,  and  the  glucocorticoid  receptor,  components  of  apparently  quite 
diverse  biological  systems,  had  some  meaningful  common  manifestation. 


-18- 


Resonant  Mode-Induced  Stimulation  and  Desensitization 
of  a  T4-Lymphocyte  Receptor 
by  a  Viral  Complex  Hydrophobic  Mode 

Ubiquitin,  a  peptide  that  becomes  covalently  conjugated  to  proteins  marked  for  non-lysosomal 
proteolytic  destruction,76  has  been  found  to  be  the  reactive  species  in  the  lymphocyte  cell  outer  mem¬ 
brane  receptor.77  Known  to  regulate  the  antigen-response  competence  of  T-  and  B-cells  via  their 
maturation  rates,  ubiquitin,  first  isolated  and  sequenced  from  thymus  by  the  Goldstein  group  (see  78  for 
review)  also  plays  a  prominent  role  in  proteolytic  systems  active  in  the  increased  rate  of  degradation  of 
abnormal  proteins.  A  second  functional  aspect  of  the  7  to  9-residue  hydrophobic  mode  of  the 
ubiquitin-containing  lymphocyte  membrane  receptor  comes  through  its  spectral  ismorphism  with  the 
glucocorticoid  receptor  and  the  known  cytotoxic  and  immunosuppressive  actions  of  adrenal  cortical  hor¬ 
mones  acting  on  T4  lymphocyte  subsets.79  There  are  examples  of "  mode-resonant"  desensitization  of 
protein  receptors  on  lymphocytes  by  high  or  prolonged  levels  of  receptor  peptide  ligands  associated  with 
their  internalization  and  degradation.80  81 

The  immunological  failure,  T4-cell  depletion,82  and  high  incidence  of  neoplastic  transformation  in 
ectoderm  and  mesodermal  tissues  that  characterize  the  acquired  immunodeficiency  syndromes  (AIDS) 
suggested  to  us  the  possibility  that  the  envelope  sequences  of  the  three  responsible  retroviruses  might 
contain  7  to  9-residue  hydrophobic  modes  as  in  the  lymphocyte-ubiquitin-EGF-glucocorticoid  receptor 
isometry.  Recent  reports  of  Addisonian  signs  of  adrenal  exhaustion  and  gland  atrophy  along  with 
plasma  ACTH  elevation  in  patients  with  AIDS  made  us  wonder  if  the  viral  peptide  envelope  might 
have  activated  and  desensitized  the  receptors  through  a  pseudo-corticoid  mode  resonance.8*-85  It  is 
known  that  hydrocortisone  greatly  facilitates  the  isolation  of  retroviruses  and  their  transmission  to  fresh 
leukocytes  from  AIDS  patients.86 

31-70  HTLV-III  . . . . . 

ARV-2  . . . . . 

LAV  . .............  . . .  .  . 
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Fig.  18  The  emergence  of  7  to  9-residue  hydrophobic  mode  wavelengths  in  the  envelope  proteins  from  three 
AIDS-related  viruses.  Amino  acid  hydrophobicities  of  HTLV-III  were  used  as  the  reference  (*);  (  +  )  and  (-) 
show  higher  or  lower  kcal  for  corresponding  residues  on  the  other  two  viruses.  In  segment  31-70,  which  fol¬ 
lows  the  putative  signal  sequence,  ARV-2  and  LAV  are  practically  identical  to  HTLV-III,  and  similarity  in 
hydrophobic  mode  spectra  would  be  expected.  However,  similar  spectra  emerge  in  the  extracellular  domain 
between  231  and  270,  where  they  differ  from  HTLV-HI  and  from  one  another.  A  7  to  9-residue  wavelength 
appears  between  591  and  630  where  there  is  virtual  identity  between  HTLV-III  and  ARV-2,  which  differ 
markedly  from  LAV.  The  prevalance  of  (  +  )  and  (-)  precludes  hydrophobically  equivalent  substitutions  of 
amino  acids  as  a  basis  for  the  spectral  similarities  and  confirms  the  complexity  of  hydrophobic  mode  transfor¬ 
mation. 
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Complex  hydrophobic  transformation  showed  7  to  9-residue  spectral  modes  in  the  envelope  pro¬ 
teins  of  the  AIDS-related  viruses  HTLV-III,  ARV-2,  and  LAV*7"89  in  the  post-signal  sequence  segment 
and  two  other  regions  (Fig.  18).  This  led  to  a  conjecture  that  the  new  anti-glucocorticoid  steroids  such 
as  RU  48690  might  prevent,  retard,  or  perhaps  facilitate  the  hypothesized  hydrophobic  mode-resonant 
viral  stimulation  and  desensitization  of  T-lymphocyte,  EGF,  and  adrenal  cortical  membrane  receptors 
by  the  AIDS  viruses.  The  idea  is  currently  being  explored. 
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